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Abstract 

Background: Identifying permissible limits of intracellular parameters such as protein expression provides 
important information for examining robustness. In this study, we used the TEV protease-mediated induction of 
protein instability (TIPI) in combination with the genetic Tug-of-War (gTOW) to develop a method to measure the 
lower limit of protein level. We first tested the feasibility of this method using ADE2 as a marker and then analyzed 
some cell cycle regulators to reveal genetic interactions. 

Results: Using TIPI-gTOW, we successfully constructed a strain in which GFP-^^^^^Ade2 was expressed at the lower 
limit, just sufficient to support cellular growth under the -Ade condition by accelerating degradation by TEV protease. 
We also succeeded in constructing a strain in which the minimal level of GFP-^^^^^Cdc20 was expressed by TIPI-gTOW. 
Using this strain, we studied genetic interactions between cell cycle regulators and CDC20, and the result was highly 
consistent with the previously identified interactions. Comparison of the experimental data with predictions of a 
mathematical model revealed some interactions that were not implemented into the current model. 

Conclusions: TIPI-gTOW is useful for estimating changes in the lower limit of a protein under different conditions, 
such as different genetic backgrounds and environments. TIPI-gTOW is also useful for analyzing genetic interactions 
of essential genes whose deletion mutants cannot be obtained. 



Background 

Effective functioning of cellular systems requires optimal 
expression of individual proteins [1-3]. On the other 
hand, cellular systems are generally robust against 
changes in protein expression [4-6]. Identif)^ing permis- 
sible limits of intracellular parameters such as protein 
expression provides important information for examin- 
ing robustness [4-6]. 

We previously developed a method designated genetic 
Tug-Of-War (gTOW), by which we can measure the 
limit of overexpression of a target protein [7,8]. In 
gTOW, we clone a target gene with its native promoter 
into a 2 (i-based plasmid harboring leu2d. Under 
leucine-lacking (-Leu) conditions, the copy number of 
the plasmid reaches >100 copies because of the selection 
bias driven by leu2d. If the target protein expression 
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reaches the upper limit, i.e., essential cellular functions 
come to a halt when the level of protein exceeds a cer- 
tain limit, the gene/plasmid copy number must decrease 
to below the upper limit. In gTOW, we thus increase 
the expression of a protein to a limit by increasing the 
copy number and then determine the upper limit by 
measuring the corresponding copy number of the gene. 
Using gTOW, we previously measured the copy number 
limits for overexpression of 30 cell cycle regulatory 
genes in the budding yeast Saccharomyces cerevisiae and 
studied the robustness of the cell cycle and mechanisms 
that can improve the robustness of the cell cycle 
[7,9,10]. We recently measured the copy number limits 
of all yeast genes and showed that the yeast cellular sys- 
tem is robust against > 100-fold increase in the copy 
numbers of each of >80% genes [11]. In that analysis, we 
also identified 115 "dosage-sensitive genes" whose copy 
number limits were <10 [11]. 

While gTOW is a method to measure the upper limit 
of protein expression, determination of the lower limit is 
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still a question. The commonly used gene-deletion ex- 
periments reduce the expression of a target protein to 0. 
Gene-deletion experiments for all genes in S. cerevisiae 
have revealed that about 20% genes are essential for cel- 
lular growth under normal conditions [12]; other genes 
are essential for cellular growth under specific condi- 
tions such as specific environmental conditions or in 
combination with other gene knockouts [13,14]. Each of 
these proteins must have a minimal requirement level 
(i.e., the lower limit) to support cellular growth under 
each of these conditions. At present, there is no effective 
method to assess the lower expression limit of a target 
protein. 

To assess the lower limit, we have to gradually reduce 
the expression of a protein and then measure the lower 
limit by some methods. To reduce the expression of a 
gene/protein, we need to reduce the production rate of 
the transcript, increase the degradation rate of the tran- 
script, decrease the production rate of the protein, or in- 
crease the degradation rate of the protein. To identify 
the lower limit of expression, we would need to specific- 
ally decrease/increase the production/degradation rate of 
the target and evaluate the limit that would halt cellular 
function. Experimental methods that could potentially 
be used for this purpose in S, cerevisiae are as follows: 
(1) decreasing the expression rate of the transcript using 
a repressible promoter [15], (2) decreasing the transcrip- 
tion rate or increasing the degradation rate using RNA 
interference (RNAi) [16], (3) increasing the degradation 
of mRNA by the decreased abundance by mRNA per- 
turbation (DAmP) method [17,18], (4) increasing the 
rate of degradation of the transcript by recruiting the 
RNA-degrading enzyme to the 3' region of the tran- 
script [19], (5) increasing the rate of degradation of the 
protein using inducible degrons [20,21], (6) increasing 
the rate of degradation of the protein using the TEV 
protease-mediated induction of protein instability (TIPI) 
method [22,23]. 

In this study, we used TIPI in combination with 
gTOW to develop a method to measure the lower limit 
of protein expression. In TIPI, a target protein is 
expressed with a TEV protease-induced degron, an 
N-terminal sequence containing a cleavage site for the 
site-specific TEV protease. Cleavage of the protein at 
this site by the TEV protease triggers the exposure of 
the "N-degron" [24], which triggers the rapid degrad- 
ation of the target protein through the N-end rule path- 
way. We chose TIPI because, theoretically, we can 
increase the degradation rate of the target protein by in- 
creasing the expression of TEV protease using the 
gTOW scheme with a plasmid encoding TEV. We can 
then evaluate the lower limit of the target protein ex- 
pression by measuring the copy number of the TEV 
plasmid. We first tested the feasibility of this method 



using ADE2 as a marker and then analyzed some cell 
cycle regulators to reveal genetic interactions. 

Results and discussion 

Establishment of TIPI-gTOW using ADE2 as a marker 

In the gTOW method, we cloned the target gene into a 
plasmid and increased the copy number using the gen- 
etic bias for leu2d [7,8]. In this study, by increasing the 
copy number of the gene encoding TEV protease using 
the gTOW scheme, we attempted to increase the deg- 
radation rate of a target protein containing a TEV 
protease-induced degron (Figure 1). If the target protein 
expression is at the lower limit for any essential cellular 
function, the degradation rate must be less than that re- 
quired to maintain the lower limit, which will restrict 
the upper limit copy number of the gene encoding TEV. 
We designated this scheme TIPI-gTOW. 

To ensure the success of this approach, we need to ad- 
just the expression of the following: (1) the target pro- 
tein (determined by the production and degradation 
rate) to be in the range required to support the growth 
of the cell (not too high and not too low) and (2) the 
TEV protease to be in the appropriate range to induce 
degradation in order to detect the lower limit of the tar- 
get protein expression. Taxis et al used the CYCl pro- 
moter and ADHl promoter to express target proteins 
and the GALl promoter to induce TEV protease expres- 
sion [22,23]. In some cases, they could not observe ex- 
pected lethal phenotypes, probably because the induced 
degradation of the target proteins by the TEV protease 
was insufficient to reduce the protein expression to their 
lower limits [22,23]. 

We used constitutive promoters to increase the ex- 
pression of TEV protease in accordance with the in- 
crease in gene copy number. We first used the S, 
cerevisiae TEFl promoter to express TEV protease, but 
the strong expression of TEV protease from the TEFl 
promoter on the high-copy gTOW plasmid in itself 
caused a cellular growth defect (data not shown). We 
currently do no know what causes this toxicity. We then 
used the CUPl promoter. In the culture conditions we 
used (0.25 (iM Cu^""), the CUPl mRNA level was about 
10% of the TEFl mRNA level (data not shown). Expres- 
sion of the TEV protease from the CUPl promoter did 
not show any cellular growth defect, even when the gene 
copy number exceeded 100 (data not shown). We used 
an efficient version of the pl4-TEV fusion protein 
(pl4*-TEV"') [22,23]; here, we designated this "pTEV." 
We also designated the gTOW plasmid encoding the 
pTEV^ expressed from the CUPl promoter as the "TEV 
plasmid." 

To test the feasibility of TIPI-gTOW, we selected 
ADE2 as a target. ADE2 encode phosphoribosylaminoi- 
midazole carboxylase, an enzyme involved in the adenine 
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Figure 1 Scheme of TIPI-gTOW. We first constructed a strain in which the chromosomal target gene was replaced by a GFP-^^^^^target 

construct. We next introduced the TEV plasmid, a plasmid for gTOW that encodes pTEV^ expressed from the CUPl promoter. According to the 
TIPI procedure, cleavage and rapid degradation of the GFP-^°^^^target is induced by pTEV^ Using the gTOW procedure, in which the copy 
number of the TEV plasmid exceeds 100 under the -Leu condition, we can increase the expression of pTEV^ which accelerates the degradation 
of the GFP-^^^^ ^target, reducing the level of the GFP-^^^^ ^target. It is thus expected that the upper limit copy number of the TEV plasmid would 
inversely correlate with the lower limit of the GFP-^^^^ ^target. The tug-of-war between the bias to increase the copy number of Ieu2d and the 
bias to decrease the copy number of pTEV^ gene determines the plasmid copy number in the cell under the -Leu condition. It is thus possible 
to indirectly estimate the lower limit of the GFP-^^^^ ^target by measuring the copy number of the TEV plasmid. 



synthesis. Reduction of Ade2 protein can be indirectly 
monitored by assessing cellular growth defect and by ac- 
cumulation of red pigment in the absence of adenine 
[25]. We integrated ADE2 with the TEV-induced degron 
expressed from the CDC19 promoter into the chromo- 
somal ADE2 locus (Figure 1). We chose the TDegF 
degron because it is one of the strongest degrons 
[22,23]. We then introduced the TEV plasmid into the 
cell. The results of TIPI-gTOW experiments of ADE2 
are shown in Figure 2. Cellular growth deficiency and 
red colony formation were observed for cells containing 
the TEV plasmid in the -Leu-Ade condition (Figure 2A). 
We note that colonies are expected to be redder on 
SC-Ura-Leu plates than the ones on SC-Ura plates, 
but they were not (Figure 2A). We currently do not 
know the reason. Unknown interaction between the 
leucine deficiency and the red color formation might 
exist. We next measured the copy numbers of the 
TEV plasmid and found that the copy number limit 
of the TEV plasmid decreased significantly under 
the -Leu-Ade condition (Figure 2B). We thus en- 
sured that the Ade2 protein was expressed at the 
lower limit and indirectly estimated the lower limit 
from the TEV plasmid copy number. We note that 
the copy numbers of the vector were significantly dif- 
ferent in between -Ura-Leu and -Ura-Leu-Ade condi- 
tions (p < 0.01, Students ^-test). The replication and/ 
or partition of the 2 (i-based plasmid might be affected by 
the growth conditions. 



TIPI-gTOW experiments of cell cycle regulators 

We next selected 3 essential cell cycle regulators as tar- 
get genes to test if TIPI-gTOW was a feasible approach 
to assess their lower limits. We constructed the GFP- 
TDegF constructs for CDC15, CDC20, and CDC28 and 
integrated these into their chromosomal loci. At first, we 
used the full-length CDC19 promoter (729 bp) to ex- 
press the cell cycle proteins but could not obtain a 
CDC20 construct, probably because strong expression of 
CDC20 is toxic [26]. We then constructed a series of de- 
letions of the CDC19 promoter to reduce its expression 
(Additional file 1: Figure SI). The GFP expression from 
promoters with a length of 500 bp or less was markedly 
lower than that from promoters exceeding 500 bp 
(Additional file 1: Figures SIB and SIC). This reduction 
might be due to the deletion of a Tye7 binding site located 
at -523 bp to -516 bp (Additional file 1: Figure SIA). We 
thus used the full-length CDC19 promoter (729 bp), the 
CDC19-600 promoter (600 bp), and the CDC19-soo pro- 
moter (500 bp) to construct the GFP -TDegF target. 

The results of the constructions are summarized in 
Additional file 1: Table SI. GFP-^^^^^CDCi5 constructs 
with all 3 promoters were obtained. We observed reduc- 
tion in the copy number of the TEV plasmid with the 
CDC 19 -500 promoter construct, but this was not signifi- 
cant (Figure 3A, p > 0.01, Students ^-test). The lower 
limit of Cdcl5 may be lower than that produced by 
TIPI-gTOW. As described above, we could not obtain a 
GFP-^^'^^CDC20 construct with the full-length CDC19 
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Figure 2 Evaluation of the lower limit of Ade2 using TIPI-gTOW. (A) Growth of the cells harboring GFP-^^''^^ADE2 (YSMOOl) with the vector 
and the TEV plasmid on SC plates without the indicated amino acids. Six independent colonies of each strain were tested. (B) Copy numbers of 
the plasmids of GFP-^^^^'' ADE2 (YSMOOl) in SC medium without the indicated amino acids. Four independent measurements were performed, 
and the average is shown. The error bar indicates the standard deviation. p < 0.01 for Student's t-test. 



promoter or the CDC19_6oo promoter probably because 
the expression was too high. However, we obtained a 
GFP-^^'^^CDC20 integrated strain with the CDC19_soo 
promoter. With this strain, we succeeded to detect sig- 
nificant restriction of the copy numbers of the TEV plas- 
mid under both -Ura and -Ura-Leu conditions 
(Figure 3B, p < 0.01, Student s t-test). We obtained GFP- 
™^^^CDC28 integrated strains with all the 3 promoters 
and observed significant restrictions of the copy number 
of the TEV plasmid with the CDC19_6oo promoter under 
the - Ura - Leu condition (Figure 3C, < 0.01, Students 
^-test). We did obtain a GFP-^^^^^CDC28 integrated 
strain with the CDC19_soo promoter but could not ob- 
tain any transformants with the TEV plasmid, probably 
because the expression of Cdc28 from the CDC19-soo 
promoter was already close to the lower limit of Cdc28, 
even in the absence of TEV protease. 

Reduction in the protein expression of GFP-^^^^^Cdc28 
by TEV protease was confirmed by western blotting with 
anti-GFP antibody (Figure 3D). The GFP-^^^^^Cdc28 pro- 
tein expressed from the CDC 19 ^500 promoter seemed to 
be too low to be detected by western blotting (Figure 3D). 
Simultaneously, we could not detect any GFP-^^^^^Cdc20 
protein expressed from the CDC19-S00 promoter (data 
not shown). 



We further analyzed the TEV protease dependent 
cleavage and reduction of GFP-^^^^^Cdc28 by western 
blotting with anti-Cdc28 antibody (Figure 3E). The 
protein levels of GFP-™^^^Cdc28 in YSM005 {P_CDC- 
19^6oo-GFP™^^^CDC28) were highly reduced in the 
presence of the TEV plasmid compared with the vector 
control, and the level was lower under the -Ura-Leu 
condition, where the expression of pTEV^ increased, 
than the -Ura condition. Cleaved products (^^^^Cdc28) 
were also observed in the presence of the TEV plasmid. 
The summed levels of the GFP-^^^^^Cdc28 and 
*^^^^Cdc28 were comparable with the Cdc28 levels 
expressed in BY4741 (wild type). These results con- 
firmed that the target protein level was actually reduced 
using TIPI-gTOW. We also tired to detect Cdc20 by 
western blotting using anti-Cdc20 antibodies (sc-6730 
and sc-6731, Santa Cruz biotechnology). However, we 
could not detect any signal, probably because the ex- 
pression level of Cdc20 is too low to detect with these 
antibodies. 

Estimation of changes in the lower limit of CDC20 in 
gene-deletion strains 

As described above, we identified conditions under 
which the GFP-^^^^^Cdc20 expression was reduced to 
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Figure 3 TIPI-gTOW experiments of the cell cycle regulators. Copy numbers of the plasmids in the strain; (A) P_CDC19_5oo-GFP-'^'^''^^CDC15 
(YSM002), (B) P_CDC19.5oo-GFP-^'^'^^CDC20 (YSM003), and (C) P_CDa9.eoo-GFP-^'^'^^CDC28 (YSM005) in SC medium without the indicated amino 
acids. Four independent measurements were performed, and the average is shown. The error bar indicates the standard deviation. p < 0.01 for 
Student's t-test. (D) Western blotting for GFP-™'^3^Cdc28 protein in strains YSM005 {P_CDC19.eoo) and YSM004 {P_CDC19.5oo) expressed in SC 
medium without the indicated amino acids, as detected with anti-GFP antibody. (E) Western blotting for GFP-™^^^Cdc28 protein and pTEV^ in 
strains BY4741 (wild type) and YSM005 {P_CDC1 9.eoo-GFP-'^'^'"^^CDC28) expressed in SC medium without the indicated amino acids, as detected 
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the lower limit. We next tried to estimate the lower 
limits of GFP-™^^Cdc20 in gene-deletion strains using 
TIPI-gTOW. We integrated the P_CDC19_soo-GFP- 
™^^^Cdc20 construct into the chromosomal CDC20 
locus of each of 23 cell cycle regulator deletion strains 
and measured the TEV plasmid copy number. The copy 
number data is shown in Table 1. We first noticed that 
the copy numbers of the vector itself varied among dele- 
tion strains. This may be because of differences in the 
growth rate or the efficiency of replication and partition- 
ing of the plasmid among the deletion strains. 

In addition, we observed significant correlations in the 
copy numbers of the vector and the TEV plasmid under 



both -Ura (r = 0.61) and -Ura-Leu (r = 7.0) conditions 
(Additional file 1: Figure S2). We thus considered that it 
was not appropriate to simply estimate the lower limit of 
GFP-™^^^Cdc20 from the copy number of the TEV plas- 
mid in each deletion strain. Therefore, we calculated the 
ratio of the copy number of the TEV plasmid to that of 
the vector in each deletion strain and compared the ratio 
with that of the wild type to estimate the relative change 
in TEV plasmid copy number (see Methods for details). 

Figure 4A and Table 1 show the relative copy number 
change in each deletion strain. Some strains, such as 
clnlA, siclA, and cdhlA, had negative relative copy 
number change values, which indicated that the copy 
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Table 1 Copy numbers of the vector and the TEV plasmid in CDC20 TIPI-gTOW in the gene-deletion strains 



Vector copy number* 
Average SD 



TEV plasmid copy number* 



Average 



SD 



TEV/ 

vector ratio 



Relative change 
in copy number 



Wild type 


149.7 


43.4 


59.8 


19.9 


0.40 


- 


bck2A 


44.3 


3.4 


25.3 


0.8 


0.57 


0.43 


bub2A 


71.4 


8.4 


48.6 


3.5 


0.68 


0.70 


cdhlA 


66.7 


5.7 


19.6 


2.2 


0.29 


-0.27 


cin8A 


53.7 


20.0 


37.1 


5.0 


0.69 


0.73 


clblA 


57.6 


5.2 


19.2 


2.2 


0.33 


-0.16 


clbIA 


56.3 


4.3 


20.2 


7.8 


0.36 


-0.10 


clbBA 


77.5 


15.1 


27.7 


3.6 


0.36 


-0.11 


clb4A 


59.4 


4.5 


19.3 


2.1 


0.32 


-0.19 


clbSA 


98.0 


21.4 


33.6 


6.5 


0.34 


-0.14 


clb6A 


76.1 


8.6 


31.6 


2.7 


0.42 


0.04 


clnlA 


80.1 


24.3 


20.7 


0.6 


0.26 


-0.35 


cln2A 


82.7 


15.8 


42.9 


3.8 


0.52 


0.30 


clnSA 


47.0 


10.9 


22.0 


10.8 


0.47 


0.17 


ItelA 


84.8 


20.9 


32.0 


16.0 


0.38 


-0.06 


mad2A 


47.4 


6.0 


35.2 


2.9 


0.74 


0.86 


mbplA 


21.6 


11.4 


15.0 


2.1 


0.69 


0.73 


mih lA 


64.4 


12.1 


29.3 


3.7 


0.46 


0.14 


Sid A 


71.8 


20.1 


20.2 


11.8 


0.28 


-0.29 


swelA 


124.4 


13.0 


42.2 


5.3 


0.34 


-0.15 


swi4A 


76.1 


15.1 


28.4 


6.6 


0.37 


-0.07 


swiSA 


94.0 


23.8 


37.3 


12.6 


0.40 


-0.01 


swi6A 


64.0 


8.7 


35.2 


7.2 


0.55 


0.38 


whiSA 


69.6 


15.8 


33.8 


10.9 


0.49 


0.22 



^Average and SD of four independent experimental trials are shown. 

numbers of the TEV plasmids in these strains were 
lower than those in the wild type, indicating that the 
lower limits of GFP-^^^^^Cdc20 increased in these 
strains. This implies that the deleted genes are involved 
in the production or activation of Cdc20 or have func- 
tions overlapping those of Cdc20 (Additional file 1: 
Figure S3). In these cases, the genes could be synthetic- 
ally lethal with cdc20 mutants or multicopy suppressors 
of cdc20 mutants. 

Other strains, such as mad2A, cinSA, and mbplA, 
showed positive relative copy number change values, in- 
dicating that the copy numbers of the TEV plasmids 
were higher than those in the wild type; thus, the lower 
limits of GFP-^^^^^Cdc20 decreased in these strains. 
These deleted genes are therefore involved in inhibition 
of Cdc20, e.g., by inhibition of Cdc20 production or 
Cdc20 activity or degradation of Cdc20 (Additional file 
1: Figure S4). These genes could be dosage suppressors 
of the toxicity of Cdc20 overexpression. 

Previous studies, including systematic analyses, identi- 
fied genes showing genetic interactions with CDC20 



[13,27-30]; these were listed in the Saccharomyces 
Genome Database (yeastgenome.org). From these inter- 
actions, we summarized relevant genetic interactions 
studied here in Table 2 and Figure 4A. Importantly, the 
relative copy number of the TEV plasmid decreased in 
deletion strains involving genes showing synthetic lethal- 
ity or dosage rescue of the cdc20 mutant, whereas the 
relative copy number of the TEV plasmid increased in 
deletion strains involving genes showing synthetic res- 
cue, phenotypic suppression, and dosage lethality of the 
cdc20 mutant. This result indicated that our TIPI- 
gTOW approach gave information regarding genetic 
interactions, which was highly consistent with the previ- 
ously identified interactions. Moreover, the advantage of 
our approach is that we could isolate genetic interac- 
tions in both directions, namely negative interactions 
such as synthetic lethal and dosage rescue and positive 
interactions such as synthetic rescue, phenotypic sup- 
pression, and dosage lethality in a single experiment. 

Among genetic interactions obtained in this study, the 
difference between CLNl and CLN2 was surprising 
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Figure 4 Evaluation of the lower limits of Cdc20 protein on deletion of cell cycle regulators. (A) Relative copy number changes in tine 
TEV plasmids in gene-deletion strains with GFP-^^^^^CDC20. Genes with previously identified genetic interactions with cclc20 mutation were boxed 
(summarized in Table 2). The original data used for the graph is shown in Table 1 . (B) Relative copy number changes in TEV plasmids in gene- 
deletion strains with GFP-^^^'^^CDC20 selected from A and simulation results. The degradation rate of kd20 is gradually increased until the cell cycle 
simulation indicates "dead." In the simulation, 1: Clbl and Clb2, 2: Clb5 and Clb6, 3: CInl and Cln2, and 4: Cln3 and Bck2 are implemented as 
single genes. 5: Swi6 is implemented into the component of SBF (together with Swi4) and MBF (together with Mbpl). The original data used for 
the graph is shown in Table 1 and Additional file 1: Table S2. 



Table 2 Known genetic interactions with CDC20 mutation listed in the Saccharomyces Genome Database^ 



Gene 


Genetic interaction with the cclc20 Mutant 


Description 


References 


BUB2 


Synthetic Rescue 


Temperature sensitivity of cdc20-l [Is) is rescued by bub2 deletion. 


[27] 


CDHl 


Synthetic Lethality 


cdc20-l its) shows synthetic lethality with cdhl deletion. 


[28] 


CIN8 


Dosage Lethality 


In the absence of CIN8, overexpression of CDC20 is lethal. 


[29] 


CLB2 


Negative Genetic 


Growth defect of cdc20_tsq368 increases with clb2 deletion. 


[13] 


LTEl 


Synthetic Lethality 


cdc20-l its) shows synthetic lethality with Itel deletion. 


[28] 


MAD2 


Phenotypic Suppression 


cdc20 blocks the premature degradation of Pdsl observed in a mod2 mutant. [30] 


SICl 


Synthetic Lethality 


cdc20-l{ts) shows synthetic lethality with sic! deletion. 


[28] 



*When synonymous genetic interactions were reported In more than one study, one representative is listed. 
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because they are known to have heavily overlapped func- 
tions [31]. Previous analysis, however, showed that CLB4 
has opposite genetic interactions with CLNl and CLN2 
[32]. Our results might indicate that these genes have 
different functions in some contexts. 

Comparison of experimental data with predictions of a 
simulation model is useful for evaluating the predictive 
ability of the mathematical model and for speculating 
about molecular mechanisms generating negative/posi- 
tive genetic interactions [7,9,10]. An integrative mathem- 
atical model of the budding yeast cell cycle was 
developed previously [33], and we have modified this 
model using gTOW data [9]. Using the "stabilization 
model" published in the reference [9], we here measured 
how much the degradation parameter of Cdc20 could be 
increased (i.e., the upper limit of kd20) in each of the 
gene-deletion models and compared the upper limit of 
kd20 with that in the wild type to obtain the relative 
change in the upper limit of kd20 (see Methods for de- 
tails). This analysis is considered to represent the process 
of TIPI-gTOW, namely indirectly assessing the lower 
limits of Cdc20 in the presence of gene deletions, by 
measuring how much degradation can increase without 
halting the cell cycle. The simulation result is shown in 
Additional file 1: Table S2. Figure 4B shows the compari- 
son of the results of TIPI-gTOW and those of the simula- 
tion. Deletion strains, such as mad2A, cdhlA, and others, 
showed good agreements between TIPI-gTOW and simu- 
lation data; however, some genes such as siclA, swiSA, and 
mbplA showed almost opposite changes, suggesting ab- 
sence of some regulatory steps in the model. 

Conclusions 

In this study, we tried to develop a method for estimat- 
ing the lower limit of protein expression by combining 
the TIPI and gTOW methods. Using TIPI-gTOW, we 
successfully constructed a strain in which GFP-^^^^^ 
Ade2 was expressed at the lower limit, just sufficient to 
support cellular growth under the -Ade condition by ac- 
celerating degradation by TEV protease. 

We also succeeded in constructing a strain in which 
the minimal level of GFP-^^^^^Cdc20 was expressed by 
TIPI-gTOW. Using this strain, we studied genetic inter- 
actions between cell cycle regulators and CDC20, We 
concluded that TIPI-gTOW is useful for estimating 
changes in the lower limit of a protein under different 
conditions, such as different genetic backgrounds and 
environments. TIPI-gTOW is useful for analyzing gen- 
etic interactions of essential genes whose deletion mu- 
tants cannot be obtained. One important characteristic 
of TIPI-gTOW is that it is useful for identifying both 
negative and positive genetic interactions, as shown in 
Figure 4. The synthetic genetic array is a potent high- 
throughput approach to identifying genetic interactions 



[34], which now allows quantitative analyses [35] and 
can also applied to essential genes [36]. TIPI-gTOW de- 
veloped in this study may be an alternative method to 
study genetic interactions in detail. 

The conventional gTOW method developed previously 
uses a gene with a native promoter and terminator as a 
unit to increase protein expression and measures the 
copy number limit to estimate the "copy number limit 
for protein overexpression" [7,8,10,11]. The advantage of 
gTOW is that we can indirectly estimate by how much a 
target protein can be overexpressed over the native level 
by simply measuring the copy number of the gTOW 
plasmid. In TIPI-gTOW, however, we first have to mod- 
ify the native copy of the target gene, which changes the 
native expression, activity, and stability of the target pro- 
tein. In addition, the aim of the procedure is to reduce 
the expression of the target protein, which makes it diffi- 
cult to quantitate the lower limit of the target protein by 
western blotting when this level is under the detection 
limit. In fact, we have failed to detect the protein levels 
of Ade2 and Cdc20 by western blotting using their spe- 
cific antibodies (data not shown). We thus believe that 
TIPI-gTOW is a method by which a cellular condition 
can be constructed with the target protein expressed at 
the lower limit still supporting cellular function. Thus, it 
is useful for identifying genetic interactions as described 
above, although it does not measure the absolute lower 
limit of the target protein. 

To use TIPI-gTOW, the optimal expression for the 
target GFP-TDegF protein must be determined so that 
the level is within the range of detection of the lower 
limit by increasing the TEV plasmid. Taxis and Knop de- 
veloped a series of TDeg constructs with different pro- 
moters and N-degrons of various strengths [23], which 
will be useful to systematically determine the optimal 
expression. 

Previous studies on haploinsufficiency and dosage 
compensation suggested that most of yeast proteins are 
expressed at least twice more than the levels required in 
diploid cells [37,38]. If we could detect the lower limit 
protein levels under the TIPI-gTOW experiments, we 
might be able to argue about this issue. We failed to de- 
tect some proteins using western blotting with specific 
antibodies as described above. Fusing the target proteins 
with more sensitive proves might be useful. 

Methods 

Yeast strains and growth conditions 

Yeast strains used in this study are listed in Table 3. 
A S. cerevisiae strain BY4741 {MATa, hisSAl, leu2A0, 
met 1 SAO, uraSAO) [39] was used as the host strain for 
the experiments. BY4741 derivatives with gene deletions 
in the cell cycle regulators were obtained from Open 
Biosystems/Thermo Scientific. YSMOOl was constructed 
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Table 3 Yeast strains used in this study 



Name 


Genotype 


Source 


BY4741 


MATa, his3A], leu2A0, metlSAO, uroSAO 


[39] 


YSMOOl 


ode2::hphMX-P_CDC19-GFP-^^'^^ 
ADE2 in BY4741 


Tiiis study 


YSM002 


cdc 1 5::hphMX-P_ CDC 1 9.500-GFP-^^'^^ 
CDC15 in BY4741 


Tliis study 


YSM003 


cdc20::hphMX-P_CDC19.5oo-GFP-^^'^^ 
CDC20 in BY4741 


Tliis study 


YSM004 


cdc28::hphMX-P_CDC19_5oo-GFP-^^'^^ 
CDC28 in BY4741 


Tiiis study 


YSM005 


cdc28::hphMX-P_CDC19.6oo-GFP-'^^'^^ 
CDC28 in BY4741 


Tiiis study 



by integrating the P_CDC19'GFP-^^'^^ADE2 fragment, 
which was amplified by PGR with the primers OHML432 
and OHML389 using pSS1002 as a template, into the 
chromosomal ADE2 locus of BY4741. Similarly, YSM002 
was constructed by integrating the P_CDC19-soo' 
QPp TDegF fj.^gj-Qgj^^^ which was amplified by PGR with the 
primers OHML504 and OHML505 using pSMOOl as 
a template, into the chromosomal CDC IS locus of 
BY4741. YSM003 was constructed by integrating the 
P_CDC19_soo'GfP-^^'^^ fragment, which was amplified 
by PGR with the primers OHML506 and OHML507 using 
pSMOOl as a template, into the chromosomal CDC20 
locus of BY4741. YSM004 was constructed by integrating 
the P_CDC19^5oo'GFP-^^'^^ fragment, amplified by PGR 
using the primers OHML508 and OHML509 and 
pSMOOl as a template, into the chromosomal CDC28 
locus of BY4741. YSM005 was constructed as for 
YSM004, except that pSM002 was used as a PGR tem- 
plate. Each P_CDC19_5oo'GfP-^^'^^CDC20 integrated 
strain harboring a gene deletion was constructed by inte- 
grating the same DNA fragment as YSM003 into the 
chromosomal CDC20 locus of each deletion strain. The 
sequences of PGR primers are listed in Additional file 1: 
Table S3. Yeast cells were cultured as described previously 
[9,40]. Yeast cells were cultured in SG media without ura- 
cil (Ura), adenine (Ade), leucine (Leu), and histidine (His) 
as indicated. 2% glucose was used as a carbon source. SG 
media were made using YNB with ammonium sulfate (MP 
biomedicals) that contains 0.25 [iM Cu^^, 

Plasmids used in this study 

Plasmids used in this study are listed in Table 4. pSS1006 
(the TEV plasmid) was constructed by inserting pl4 - 
TEV"" from pDS5 [22] into pSBI40 [7] so that pl4*-TEV"' 
is expressed under control of the CUPl promoter. pSS1002 
was constructed by connecting hph-MX4, a hygromycin 
B-resistance cassette [41], the CDC19 promoter, GFP- 
TDegF from pDS41 [22], and ADE2, followed by insertion 
into pRS423ks [11]. pSS1002, with different lengths of 
CDC19 promoters, was constructed using the PGR primers 



Table 4 Plasmids used in this study 


Name 


Relevant characteristics 


Source 


pSBI40 


2 jj on, URA3, leu 2d; plasmid 
for gTOW experiment 


[7] 


pRS423ks 


2 lA ori, HISS 


[11] 


pDS5 


pl4*-TE\A 


[22] 


pDS41 


GFP-TDegF-SFSb 


[22] 


pSS1002 


hph-MX4, P_CDCI9-GFP-'^^'^^ 
ADE2 in pRS423ks 


Tliis study 


pSS1006 


P_CUPI-PI4*-TE\A in pSBI40 


Tliis study 


pSMOOl 


P_CDCI9.5oo instead of P_CDC 19 
in pSS1002 


Tliis study 


pSM002 


P_CDCI9-6oo instead of P_CDC 19 
in pSS1002 


Tliis study 



OHML541-547 and OHML540; the sequences are listed in 
Additional file 1: Table S3. 



Copy number determination of plasmids 

The copy number of a plasmid in a cell was measured as 
described previously [7] with some modifications. Gells 
of each yeast strain with pSS1006 (the TEV plasmid) or 
pSBI40 (the empty vector plasmid) were grown for 50 h 
in 200 [iL SG medium without relevant amino acids. 
The cells were corrected and suspended into 50 [iL of 
Zymolyase solution (10 mM Na-phosphate [pH 7.5], 
1.2 M sorbitol, and 2.5 mg/ml Zymolyase 100 T (Nacalai 
tesque)). The cell suspension was incubated for 15 min 
at 37°G, and then incubated for 10 min at 100°G. After 
removing the cell debris by centrifugation, the supernatant 
was used as the total DNA solution. Two real-time PGRs 
were performed using 2 \iL of the DNA solution as tem- 
plates with LEU2 and LEU3 primer sets (primer sequences 
are listed in Additional file 1: Table S3), using LightGycler® 
480 SYBR Green I Master with the LightGycler" 480 
system (Roche). The LEU2 and LEU3 primer sets were 
used to quantify the LEU2 DNA from the plasmid 
and the LEU3 DNA from the genome, respectively. 
The plasmid copy number per haploid genome was 
estimated by comparing the relative amount of the LEU2 

DNA and the LEU3 DNA as follows: Plasmid copy 
number = 2^'^^ '''^'"'^^-"'^'^ , where Gp_LEU2 and Gp_LEU3 
are the crossing points of PGR amplifications using the 
LEU2 and LEU3 primer sets, respectively. 

The relative copy number change of the TEV 
plasmid was calculated as follows: Relative copy num- 
ber change = {ACt/Cy- wtCT/wtCyj/iwtCT/wtCv), where 
ACt is the average copy number of the TEV plasmid in 
each knockout strain, ACy is the average copy number 
of the empty vector plasmid in each knockout strain, 
wtCr is the average copy number of the TEV plasmid 
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in the wild-type strain, and wtCy is the average copy 
number of the empty vector plasmid in the wild-type 
strain. 

Western blotting 

Cells were cultivated in SC medium, and the proteins 
were extracted as described previously [7] . The proteins 
were separated using NuPAGE® 4-12% Bis-Tris gels, 
and transferred them onto a PVDF membrane using 
iBlot® gel transfer system (Life Technologies), then the 
target proteins were detected using their specific anti- 
bodies. For the detection of GPP, TEV protease, and 
Cdc28, antibodies 1814460 (Roche), PAB19931 
(Abnova), and sc-6709 (Santa Cruz Biotechnology) were 
used, respectively. For the second antibodies, Histofine® 
simple stain MAX PO (MULTI) and Histofine® simple 
stain MAX PO (G)(Nichirei) were used. Coomassie® 
G-250 (SimplyBlue™ SafeStain, Life Technologies) was 
used to stain loading controls. 

Computation 

Numerical simulations were performed using Matlab 
version 7.3.0 with the "stabilization model" described 
previously [9]. The rate of degradation of kd20 was grad- 
ually increased until the cell cycle simulation yielded the 
result "dead." The relative degradation rate change was 
calculated as follows: Relative degradation rate change = 
{AULkd2o-'^tULkd2o)l'^tULM2o^ where is the 

upper limit of kd20 in a deletion mutant model and 
wtUL]^d2o is that in the wild- type model. 

Additional file 
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